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Abstract. We present the results from a 500 ks Chandra observation of the Seyfert 1 galaxy NGC 5548. We detect broadened 
(full width half maximum = 8000 km s"^) emission lines of O VII and C VI in the spectra, similar to those observed in the 
optical and UV bands. The source was continuously variable, with a 30 % increase in luminosity in the second half of the 
observation. The gradual increase in luminosity occurred over a timescale of ~ 300 ks. No variability in the warm absorber 
was detected between the spectra from the first 170 ks and the second part of the observation. The longer wavelength range 
of the LETGS resulted in the detection of absorption lines from a broad range of ions, in particular of C, N, O, Ne, Mg, Si, S 
and Fe. The velocity structure of the X-ray absorber is consistent with the velocity structure measured simultaneously in the 
ultraviolet spectra. We find that the highest velocity outflow component, at —1040 km s~^, becomes increasingly important 
for higher ionization parameters. This velocity component spans at least three orders of magnitude in ionization parameter, 
producing both highly ionized X-ray absorption lines (Mg XII, Si XIV) as well as UV absorption lines. A similar conclusion is 
very probable for the other four velocity components. 

Based upon our observations, we argue that the warm absorber probably does not manifest itself in the form of photoionized 
clumps in pressure equilibrium with a surrounding wind. Instead, a model with a continuous distribution of column density 
versus ionization parameter gives an excellent fit to our data. From the shape of this distribution and the assumption that the 
mass loss through the wind should be smaller than the accretion rate onto the black hole, we derive upper limits to the solid 
angle as small as 10~^ sr. From this we argue that the outflow occurs in density-stratified streamers. The density stratification 
across the stream then produces the wide range of ionization parameter observed in this source. We determine an upper limit of 
0.3 M0 yr"^ for the mass loss from the galaxy due to the observed outflows. 

Key words. AGN: Seyfert 1 galaxy - X-ray: spectroscopy - quasar: individual: NGC 5548 



1. Introduction 

Over half of all Seyfert 1 galaxies exhibit signatures of pho- 
toionized outflowing gas in their X-ray and UV spectra. 
Studying these outflows is important for a better understand- 
ing of the enrichment of the Inter Galactic Medium (IGM) as 
well as the physics of accretion of gas onto a super-massive 
black hole. 
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One of the best candidates to study these processes 
is NGC 5548. A recent XMM-Newton Reflection Grating 
Spectrometer (RGS) observation of NGC 5548 showed evi- 
dence for an X-ray warm absorber that has a nearly continu- 
ous distribution of ionization states (Steenbrugge et al. l2003t . 
There is a pronounced increase in column density as a func- 
tion of ionization parameter The average outflow velocity mea- 
sured in the X-ray band is consistent with the range of outflow 
velocities measured in the UV band, however, the spectral res- 
olution in the X-ray band is not sufficient to resolve the in- 
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dividual velocity components observed in the UV band. The 
detection of the ls-2p line of Ovi in X-rays allows for a di- 
rect comparison with the column density as determined from 
measurements of the 2s-2p lines by FUSE in the UV band 
(Brotherton et al. 2002 1. From this direct comparison, as well as 
the total column density derived from other lowly ionized ions, 
Arav et al. (2003 1 and Steenbrugge et al. ( I2003> concluded 
that there is substantially more lowly ionized gas than has been 
claimed from previous UV observations. It was concluded that 
the X-ray and UV warm absorbers are different manifestations 
of the same phenomenon. Mathur, Elvis & Wilkes ( 1995 1 al- 
ready noted this possibility before the advent of high resolution 
X-ray spectroscopy. However, their model included only one 
ionization component and this is unrealistic for the observed 
NGC 5548 high resolution spectra. 

The difference in derived column densities in the UV band 
versus X-rays can be at least partially explained if the UV ab- 
sorption lines do not cover the naiTow emission line (NEL) re- 
gion. The column densities derived from the UV lines are then 
only lower limits due to saturation effects, because these lines 
then locally absorb almost all the radiation from the continuum 
and the broad emission line (BEL) region (Arav et al. 120021 
Arav et al. l2003t . 

The discrepancy in the Ovi column density mentioned 
above is based on non-simultaneous X-ray and UV observa- 
tions and therefore it remains possible that the difference is 
due to variability in the absorber The present Chandra obser- 
vation was proposed to remedy this, by having simultaneous 
HST STIS, FUSE and Chandra observations. Unfortunately, 
due to technical problems FUSE was unable to observe at the 
scheduled time, so that no simultaneous O VI observations were 
obtained. However, the spectral resolution of Chandra allows 
us to resolve the — 1040 km s^^ component from the four other 
velocity components detected in the UV, for the strongest lines. 
This is a rather stringent test of the kinematic relation between 
both absorbers. 

The long wavelength range of the LETGS allows us to map 
out the column density as a function of ionization parame- 
ter. Using only iron to avoid abundance effects we can sample 
about three orders of magnitude in ionization parameter from 
Fe VI to Fe XXIV. Due to the high signal-to-noise ratio, we are 
sensitive to changes in ionization parameter as small as 0. 15 in 
log This is the first Seyfert 1 spectrum where the signal to 
noise is such that we are able to study in detail the absorption 
lines between 60 — 100 A. 

The details of the observations and the data reduction are 
given in Sect. 2. In Sect. 3 we present the spectral data analy- 
sis. The short and long term spectral variability of the warm ab- 
sorber are discussed in Sect. 4. In Sect. 5 we discuss our results. 
In Sect. 6 we discuss the ionization structure, the geometry and 
the mass loss through the outflow discussed in Sect. 7. The con- 
tinuum time variability is discussed in detail in Kaastra et al. 
(2004) . Limits on the spatial distribution of the X-ray emis- 
sion are given by Kaastra et al. J2003> . The data were taken 
nearly simultaneously with HST STIS observations detailed by 
Crenshaw et al. (i2003j. 



Table 1. The exposure time and the observation details for the 
Chandra and HST observations of NGC 5548. ACIS stands for 
the Advanced CCD Imaging Spectrometer and HRC for High 
Resolution Camera. The E140M grating covers the wavelength 
range between 1150 — 1730 A, the E230M grating between 
1607 -3119 A. 



Instrument 


detector 


start date 


exposure (ks) 


HETGS 


ACIS-S 


2002 Jan. 16 


154 


LETGS 


HRC-S 


2002 Jan. 18 


170 


LETGS 


HRC-S 


2002 Jan. 21 


170 


HST STIS 


E140M 


2002 Jan. 22 


7.6 


HST STIS 


E230M 


2002 Jan. 22 


2.7 


HST STIS 


E140M 


2002 Jan. 23 


7.6 


HST STIS 


E230M 


2002 Jan. 23 


2.7 



2. Observation and data reduction 

NGC 5548 was observed for a full week with the Chandra High 
Energy Transmission Grating Spectrometer (HETGS) and the 
Low Energy Transmission Grating Spectrometer (LETGS) in 
January 2002. In Tabled the instrumental setup and the expo- 
sure times are listed. 

The LETGS data were reduced as described by Kaastra et 
al. J2002at . The LETGS spans a wavelength range from 1 — 
180 A with a resolution of 0.05 A (full width half maximum, 
FWHM). The data were binned to 0.5 FWHM. However, due to 
Galactic absorption toward NGC 5548, the spectrum is heavily 
absorbed above 60 A, and insignificant above 100 A. In the 
present analysis we rebinned the data between 60 — 100 A by 
a factor of 2 to the FWHM of the instrument and ignore the 
longer wavelengths. 

The HETGS data were reduced using the standard 
CIAO software version 2.2. The HETGS spectra consist 
of a High Energy Grating (HEG) spectrum which covers 
the wavelength range of 1.5 - 13 A with a FWHM of 
0.012 A, and a Medium Energy Grating (MEG) spectrum 
which covers the wavelength range of 1.5 — 24 A with a 
FWHM of 0.023 A. The MEG and HEG data are binned 
to 0.5 FWHM. In the reduction of both HETGS datasets 
the standard ACIS contamination model was applied (see 
http://cxc.harvard.edu/caldb/about_CALDB/index.html). 

The data shortward of 1 .5 A were ignored for all the instru- 
ments, as were the data above 11.5 A for the HEG spectrum 
and above 24 A for the MEG spectrum. The higher orders are 
not subtracted from the LETGS spectrum, but are included in 
the spectral fitting, i.e. the response matrix used includes the 
higher orders up to order ±10, with the relative strengths of 
the different orders taken from ground calibration and in-flight 
calibrations of Capella. We did not analyze the higher order 
HETGS spectra, due to their low count rate. All the spectra 
were analyzed with the SPEX software (Kaastra et al. l2002bt . 
and the quoted errors are the 1 a rms uncertainties, thus Ax^ = 
2. 

Kaastra et al. (:2002al found a difference in the wave- 
length scale between the MEG and LETGS instrument of about 
170 km s^^. This wavelength difference is within the absolute 
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calibration uncertainty of the MEG instrument. The best line in 
NGC 5548 to test the wavelength calibration is the O VII for- 
bidden emission line, as it is found to be narrow. However, due 
to the noise in the MEG spectrum, the wavelength is not well 
determined, but it is consistent with the wavelength measured 
from the LETGS spectrum. We also looked for differences in 
the line centroid of the O VIII Lya and O VIII Ly/3 absorption 
line. However, as these lines are broadened, centroids are less 
accurate. These lines have consistent centroids for both instru- 
ments, and no systematic blue- or redshift is detected. 

3. Spectral data analysis 

In the online appendix C the spectrum between 1.5 A — 100 
A is shown. The X-ray spectrum is qualitatively similar to 
the spectra observed before by Chandra (Kaastra et al. 2000 
Kaastra et al. I2002al Yaq oob et al. l200Tt and XMM-Newton 
(Steenbrugge et al. '2003V A large number of strong and weak 
absorption lines cover a smooth continuum. In addition, a few 
narrow and broadened emission lines are visible. Most absorp- 
tion lines can be identified with lines already observed in these 
previous spectra or with predicted lines that, due to the in- 
creased sensitivity of the present spectra, become visible for 
the first time. 

The spectral model applied here is based on results ob- 
tained from earlier Chandra and XMM-Newton spectra of 
NGC 5548 (Kaastra et al. l2002al Steenbrugge et al. 2003 1. The 
current spectrum is described by the same spectral components 
as in the earlier observations, albeit with different parameters. 
In the next subsections we discuss the results for the differ- 
ent components as obtained from our global fit to the spectra. 
These previous observations as well as the present observation 
have a continuum that is well fitted by a power-law and mod- 
ified black body component (Kaastra & Barr .l989j modified 
by Galactic extinction and cosmological redshift. The Galactic 
Hi column density was frozen to 1.65 x 10^"* m~^ (Nandra 
et al. [T993 ) . as was the redshift of 0.01676 from the optical 
emission Unes (Crenshaw & Kraemer ll999> . The continuum 
parameters for our best fit model, namely the one with three 
slab components (Sect. l3^m model A), are listed in Table |2] 
For comparison, the continuum parameters for the earlier RGS 
observation are included (Steenbrugge et al. "2003V 

The high resolution spectra reveal a warm absorber (WA) 
modeled by combinations of slab, xabs or warm components 
(Sect. B.lt . narrow emission lines (Sect. l3.4l and l3.5> and broad- 
ened O VIII, O VII and C vi emission lines (Sect. l3.6> . 

3.1. Warm absorber models 

In our modeling of the warm absorber we use combinations 
of three different spectral models for the absorber The most 
simple is the slab model (Kaastra et al. .2002a i. which calcu- 
lates the transmission of a slab of a given column density and 
a set of ion concentrations. Both the continuum and the line 
absorption are taken into account. As in the other absorption 
models, the lines are modeled using Voigt profiles. Apart from 
the ionic column densities, the average outflow velocity v and 
the Gaussian r.m.s. velocity broadening a are free parameters 




100 1000 
Wavelength (A) 



Fig. 1. The SED used (solid line) in the present analysis and 



the analysis of the HST STIS data by Crenshaw et al. ( |2003> . 
For comparison, the SED used by Kaastra et al. ( 2002at and 
Steenbrugge et al. (2003( in earlier papers on NGC 5548 is 
plotted as the dot-dash line. The SED adopted by Dumont et al. 
( 1998 ) is plotted as the dotted fine. 

Table 2. The best fit parameters for the continuum for the ear- 
lier RGS observation (Steenbrugge et al. 2003 1 and the present 
HETGS and LETGS observations of NGC 5548. The parame- 
ters for the HETGS and LETGS are for the best fit model in- 
cluding the three slab components with the outflow velocities 
frozen to the UV values. 



Parameter 



RGS 



HETGS 



LETGS 



pi: norm" 

phlum*" 

pl:r 

mbbmorm'^ 
mbbiT (in eV) 



1.57 ±0.03 0.90 ±0.01 1.30 ±0.01 

5.7 ±0.1 3.62 ±0.04 4.02 ± 0.03 

1.77 ±0.02 1.71 ±0.01 1.88 ±0.01 

6.5 ± 0.7 < 400 3.5 ± 0.9 

97 ± 6 - 100 ± 10 



" At 1 keV in 10^^ ph s"^ keV"\ 
The luminosity as measured in the 2 — 10 keV band in 10'^'' W. 
norm = emitting area times square root of electron density in 10"^ 



of the model. The last two parameters have the same value for 
each ion. The slab is assumed to be thin in the sense that e^^ 
gives a fair description of its transmission. 

The second model (xabs) is similar to the slab model, how- 
ever now the ionic column densities are coupled through a set 
of photoionization balance calculations (Kaastra et al. l2002at . 
As a result, instead of the individual ionic concentrations we 
now have the ionization parameter ^ = Lln^^r^ and the ele- 
mental abundances as free parameters. The set of photoion- 
ization balance calculations was done using the XSTAR code 
(Kaflman & Krolik 1999 1, with the SED as shown in Fig.[T] 
The same SED was used by Crenshaw et al. J2003> in the anal- 
ysis of the simultaneous HST STIS data. As there is no evi- 
dence of a blue bump this was not included in the currently 
used SED. This is the main difference with the SED used by 
Kaastra et al. (2002a) and Steenbrugge et al. J2003I I. which is 
plotted as the dot-dash line. For the X-ray part of the spectrum 
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we used the continuum model as measured from the present 
LETGS spectrum, instead of the continuum obtained from the 
1999 LETGS spectrum. For the high energy part of the SED 
we used a power law component as measured in our Chandra 
spectrum but with an exponential cut-off at 130 keV, consistent 
with the BeppoSAX data (Nicastro et al. 2000 1. As reflection 
components are time variable and the LETGS has insufficient 
high energy sensitivity to model it properly, we did not include 
any reflection component in our SED. Between 50 and 1000 
A we used a power-law and at even longer wavelength used the 
energy index obtained by Dumont et al. (Q^998 1 of 2.5. 

The third method used is the warm model, which is a model 
for a continuous distribution of the dA'^n/d^ as a function of ^. 
At each value of ^, the transmission is calculated using the xabs 
model described above. In our implementation of this model 
^dA^n/d^ is defined at a few (in our case two) values of ^. 
At the intermediate points the logarithm of this function is de- 
termined by cubic spline interpolation on the log ^ grid. This 
determines A'h as a function of ^. Integration over ciA^n/d^ be- 
tween the lowest and the highest values of ^ (^i respectively 
^2) yields the total ionic column densities. In our model, as we 
have taken only two points, the total hydrogen column den- 
sity is therefore a power-law function of ionization. The warm 
model correctly incorporates the fact that ions are formed at a 
range of ionization parameters. Free parameters of this model 
are Ci. 6. dA^n (Ci)/dC, 6 dA^n (C2)/dC, and similar to the 
xabs model the outflow velocity v, velocity broadening a and 
the elemental abundances. 

Finally, in all our models (slab, xabs and warm) we 
changed the wavelength for the O V ls^2s^ - ls2s^2p ^Pi X-ray 
fine from 2 2.33 A (the HULLAC value) to 22.374 A (Schmidt 
et al. I2004> as measured at the University of California EBIT-I 
electron beam ion trap. This recent value was not yet avail- 
able for our earlier analysis of the RGS data (Steenbrugge et 
al. l2003 l. 

3.2. Spectral fits with warm absorber models 

3.2.1 . Warm absorber model A: column densities 
with UV velocity structure 

Our first method (dubbed model A here) allows us to further as- 
certain that the X-ray warm absorber and the UV absorber are 
the same phenomenon, a major goal for proposing this observa- 
tion. We fitted the warm absorber in the X-ray spectra with the 
outflow velocities frozen to those measured in the UV spectra. 

We implemented this by modeling the WA using three slab 
components, with the outflow velocities frozen. These outflow 
velocities were frozen to —1040 km s^^, —530 km s^^ and 
— 160 km s~^, corresponding to UV velocity components 1, 
the average of 2 — 4, and 5, respectively. Note that the veloc- 
ities of components 2 — 4 (—667, —530 and —336 km s^^, 
respectively) are too close to be separated in our X-ray spectra. 
In the fit we left the velocity broadening of each of the three 
components a free parameter. Allowing for three velocity com- 
ponents instead of one, substantially improved the fit from 
= 2916 to 2441 for 2160 degrees of freedom (only LETGS is 
quoted here). However, the absorption fines are still only par- 



Table 3. Parameters for six ions measured using model A. The 
first row gives the component number as listed in the UV. The 
second row lists the outflow velocity v, which was frozen to the 
UV values, the third row lists the velocity broadening ax-ray 
as derived from the X-ray spectra. In the fourth row the ve- 
locity broadening ctuv = FWHM/2.35, where FWHM is the 
measured and resolved Full Width Half Maximum of the UV 
absorption lines (Crenshaw et al. l2003> . All are in km s^^. For 
the —530 km s~^ component we added the velocity broaden- 
ing on component 3 and 4. We list the logarithms of the column 
densities in m~^. For comparison we also list the C IV and N V 
column densities as measured by Crenshaw et al. (|2003> . 



comp. 


1 


2-4 


5 




V 


-1040 


-530 


-160 




0"X-ray 


40 ±5 


100 ± 15 


90 ± 13 




o"uv 


94 ±8 


140 ± 15 


26 ± 6 




ion 








total 


Civ 


18.05 ± 0.05 


18.66 ± 0.02 


17.76 ± 0.03 


18.8 


Cvi 


20.2 ± 0.6 


21.4 ± 0.2 


20.8 ± 0.4 


21.5 


Nv 


18.44 ± 0.02 


19.24 ± 0.02 


18.16 ±0.03 


19.3 


Ov 


20.0 ± 0.5 


20.5 ± 0.3 


20.2 ± 0.6 


20.8 


Ovii 


21.3 ±0.5 


21.4 ± 0.3 


20.3 ± 0.7 


21.7 


VIII 


22.2 ± 0.1 


21.5 ±0.3 


21.9 ±0.3 


22.4 


Ne IX 


20.5 ± 0.9 


20.4 ± 0.6 


20.7 ± 0.8 


21.0 


Si XIII 


20.8 ± 0.6 


20.6 ± 0.6 


20.5 ± 1.1 


21.1 



tially resolved, leading to strongly correlated errors for the de- 
rived column densities of the velocity components. Only for the 
six ions with the deepest lines (C VI, O V, O VII, O VIII, Ne IX 
and Si Xlll) we measure column densities with meaningful er- 
ror ranges (Table |3}- For these six ions the velocity structure 
of the X-ray warm absorber closely resembles the UV velocity 
structure. Fig.|2]shows the line profiles as a function of outflow 
velocity for four lines of the ions listed in Table |3] The HEG 
spectrum shows most clearly substructure in the Mg XII Lya 
absorption lines (Fig.[2|in Appendix C online). 

3.2.2. Warm absorber model B: column densities 
with simple velocity structure 

Free velocity broadening 

Using model A we can only derive accurate column densi- 
ties for six ions. In order to obtain column densities for other 
ions we need to reduce the number of free parameters. To do so 
we use a model with one outflow velocity component, leaving 
the velocity broadening as well as the outflow velocity free. 
This allows us to obtain column densities, outflow velocities 
and velocity broadening for the ions with stronger lines. We 
first made a fit with a single slab component. This slab compo- 
nent has a single outflow velocity and velocity broadening for 
all the ions. After this fit is obtained, we freeze all the parame- 
ters in this slab component. 

A second slab component is then added containing only 
a single ion, leaving the column density, outflow velocity and 
velocity broadening as free parameters. The column density for 
this particular ion is frozen to zero in the first slab component. 
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Table 4. The best fit column densities using model B for the ions for which the outflow velocity and velocity broadening are 
well constrained. The outflow velocity v and the r.m.s. velocity broadening a are listed as well as the ionization parameters f and 
U for which the ion has its maximum column density. The column densities are from the fit including the broad emission lines 
(Sect.llet. 



ion 


log iVion 


— V 




1,-,,, t 


log u 




(m ) 


{KIR S ) 


{KIR S ) 


(iU w 




L, V 


01 nine 
Zi.Z iiz U.J 


^■^n -1- on 
jjV ± yu 


LoJ =h /U 


n 1 n 

U. lU 


1 '^n 

— 1 . jU 


\^ VI 


ZL.D ± U. I 


A on _L 
4oU ± OU 


1 on _!_ 
loU d= jj 


1 nn 
I.UU 


n An 
— U.OU 


JN VI 


ZL.Z rh U.Z 


^^n _L 1 in 
jjU ± IjU 


onn _!_ ion 
ZUU dz IZU 


n ^n 
U.jU 


1 in 
— I.IU 


N VII 


Zl.j ± \J.Z 


^on 1 1 1 n 
jZU ± 1 iU 


^Qn -1- 1 on 
jy\J zt iZU 


i . J J 


n on 

—U.ZU 


U V 


on Pino 


^Qn _L 1 on 
jjU ± IZU 


oon _!_ 1 nn 
ZoU dz lUU 


n on 
— U.ZU 


1 on 


U VI 


on A -1- n 
ZU.o ± U.Z 


^cn -1- 1 ^n 
jJSU ± IjU 


1 An 1 1 1 n 
lOU zt i lu 


n ^n 
U. jU 


1 ^n 


U Vll 


00 1 Q -L n m 
ZZAo ± U.Uj 


^on -1- An 
jVU ± OU 


1 '^n -1- on 
1 jU zt ZU 


n Q'^ 

U.Vj 


n A^ 


VIII 


oo _L n n/1 
ZZ.j3 zh U.U4 


^/in _L *^n 
j4U ± OU 


1 ^n _!_ on 
I jU zt ZU 


1 A^ 
I.OJ 


n n< 
U.Uj 


Ne IX 


01 1 n ^ 

Zl.Z it U.J 


AAn -1- 1 "yn 
OOU ± 1 /U 


in -1- 1 1 n 
ZlU zt 1 lU 


1 ^n 


n 1 n 

— U. lU 


Ne X 


o 1 n 1 n o 
zi.y ± U.Z 


QQn _L 1 /in 
ojU ± 14U 


onn _!_ An 
zyu zt oU 


o n^ 
Z.Uj 


U.4j 


Mg IX 


21.0 ± 0.3 


560 ± 70 


55 zb 30 


1.00 


—0.60 


Mgxii 


21.3 ± 0.1 


680 ± 120 


430 zb 120 


2.35 


0.75 


Si VIII 


20.2 ± 0.2 


620 ± 210 


330 zb 190 


0.40 


-1.20 


Six 


20.2 ± 0.2 


790 ± 150 


340 zb 150 


1.05 


-0.55 


Si XI 


20.0 ± 0.2 


810 ± 180 


350 zb 150 


1.50 


-0.10 


Si XIII 


21.0 ± 0.1 


660 ± 280 


80zb40 


2.20 


0.60 


Si XIV 


21.0 ± 0.3 


880 ± 120 


60zb30 


2.60 


1.00 


Sxi 


20.0 ± 0.3 


560 ± 200 


310 zb 170 


1.20 


-0.30 


Sxii 


20.2 ± 0.2 


620 ± 300 


450 zb 200 


1.45 


-0.15 


Fexvil 


20.3 ± 0.4 


740 ± 110 


50zb30 


1.80 


+0.20 




Fig. 2. The MEG (open circles) and LETGS (filled squares) line 
profiles for the deepest line of O VIII, O VII, O V and C VI. The 
five outflow velocities measured in the UV are indicated by the 
dotted line. No correction was made for possible blending. 



and the spectra are fit again in order to fine tune the column 
density, outflow velocity and velocity broadening for each ion. 
Each ion was cycled through the second slab component in this 
manner, and the results for those ions which have well deter- 
mined outflow velocity and velocity broadening are listed in 
Table 13 In this Table we also list the ionization parameter for 
which the column density of the particular ion peaks. We list 
^ = L/rigr^ expressed in 10^^ W m, where L is the energy lu- 
minosity, fie is the electron density, r the distance. We also list 



□ c 

o N 



A Ne 
Mg 



0.5 1 1.5 2 2.5 

Log f (10"^ W m) 

Fig. 3. The measured outflow velocity versus the logarithm 
of the ionization parameter. The dotted line indicates the 
—530 km s^^ outflow velocity, the dominant component in 
the UV absorber. A trend toward higher outflow velocities for 
higher ionized ionization parameters is observed. 



U = (5tot/(47rr^cne), with Qtot/(47rr^) the total photon flux 
for photons with an energy greater than 13.6 eV. This Table 
indicates a weak correlation between outflow velocity and the 
ionization of the ion, which is shown in Fig.lU As a side effect, 
leaving the outflow velocity free potentially allows to correct 
for any remaining wavelength inaccuracies. 
Fixed velocity broadening 
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As it is not possible to determine the outflow velocity and 
broadening accurately for all ions, we also fitted the spectra 
using a fixed velocity broadening to further reduce the num- 
ber of free parameters. The velocity broadening was fixed to 
140 km s^^ (average a for Ovil and Ovill) and 70 km s^^ 
(average a for Si XIII and Si Xiv). Those two values were cho- 
sen to study the influence of saturation and velocity broadening 
on the determined column density. In this fit the outflow veloc- 
ity was fixed to either the derived value (see Table|3 or a value 
between —400 and —800 km s^^, depending on the ionization 
parameter of the ion. This method allows us to obtain an accu- 
rate column density for the less abundant ions. The results of 
this fit are listed in the appendix A as TablellSI As in Table|4]we 
list the ionization parameter for which the column density of 
the particular ion peaks. For several ions we find a higher col- 
umn density and larger velocity broadening in Table |4] which 
can be explained by the fact that a velocity broadening of 140 
km s~^ does not include the full blend due to the five different 
velocity components. 

3.2.3. Warm absorber model C: multiple separate 
ionization components 

In our previous models A and B we determined ionic column 
densities independently from the photoionization balance. A 
disadvantage of this method is that the column densities of the 
many ions with small abundances are poorly constrained, while 
the combined effect of these ions may still be noticeable in the 
spectrum. Therefore we have tried spectral fits with three xahs 
components. 

The fit with three xahs components is not statistically ac- 
ceptable. As in previous papers the iron abundance is mea- 
sured to be overabundant for lower ionization ions (Kaastra et 
al. l2002al Blustin et al. 12001 Steenbrugge et al. l20U3l Netzer 
et al. I2003I I. A possible explanation was given by Netzer et 
al. J2003I I. Namely, the ionization balance for M-shell iron in 
low temperature, photoionized plasmas depends strongly on 
the dielectronic recombination rates. These rates are sometimes 
poorly known due to resonance effects at low energies, and are 
thus not correctly incorporated in the photoionization balance 
codes that are used to predict the ionic column densities. In 
recent articles Kraemer, Ferland & Gabel (2003 1 and Netzer 
J2003> estimate that the ionization parameter for the iron ions 
that produce the UTA can change by as much as a factor of two 
or more. An increase of the effective ionization parameter ^ by 
a factor ^ 2 for Fe VII - Fe XIII indeed would bring these ions 
more in line with the trend observed for the other elements. 

Fitting iron simultaneously with the other elements will 
bias the ionization parameter measured. A solution is to de- 
couple iron from the other elements. Since however we need 
at least three xahs components, this would imply doubling of 
the number of free parameters to six xahs models. Another op- 
tion is to fit instead iron separately with a slah model. The re- 
sults of this fit are listed in Table |5l However, as iron has the 
largest ionization range and is the only element that samples 
both the highest ionized and lowest ionized absorber, both op- 
tions severely hmit the ionization range we are able to detect 



Table 5. The best fit results for a model with three xahs com- 
ponents fitting all elements but iron, which was fit separately 
with a slah model (model C). To fit the full blend we froze the 
velocity broadening to 200 km s^^ and the outflow velocity to 
- 530kms-i. 



Comp. 


A 


B 


C 


LETGS: 








log Nh (m-^) 


< 25.9 


< 25.8 


24.56±0.09 


log^ (10~^ Wm) 


< 5 


2.1±0.6 


0.8±0.1 


LETGS, HEG and MEG: 








log Nh (m^^) 


< 28 


25.43±0.06 


24.55±0.08 


log^ClO"'' Wm) 


< 5 


2.17±0.05 


0.80±0.09 



unbiased. The error bars derived for the fit without iron are 
much larger; and the ionization parameters measured should 
not be compared with those listed in our earlier RGS analysis, 
as we included iron in the fit. Forcing a large velocity broad- 
ening to fit all the outflow velocity components, and fitting the 
iron ions separately with a slah model, we find a decent fit for 
three ionization components (x^ = 3435 for 2758 degrees of 
freedom). Adding more xahs components leads to highly cor- 
related errors and fitting results. We list the results of this fit in 
TablelU 

Considering the error bars, due to the exclusion of iron in 
the fit, the model with different ionization parameters is only 
constrained for two ionization components (B and C) using the 
full observation. In order to further reduce the number of free 
parameters, and thus to obtain a better constrained fit as well as 
to test whether the ionization parameter distribution is contin- 
uous, we next test a model with a continuous ionization struc- 
ture. 

3.2.4. Warm absorber model D: continuous 
ionization structure 

Steenbrugge et al. (2003) found a power law like distribution 
of deduced hydrogen column density versus ionization param- 
eter. In order to test whether this apparent power law distribu- 
tion can be produced indeed by an underlying power law, we 
have fitted the spectrum directly to such a distribution using the 
warm model (Sect. lTTl . Similar to model C, we get a poor fit 
if iron is included. Modeling all iron ions using a slah model 
and all remaining elements with the warm model, we obtain 
a good fit to the data, namely = 3305 for 2754 degrees of 
freedom, assuming solar abundances. Table|6llists the best fit 
parameters using the warm model. In this fit the best fit outflow 
velocity is —530 km s^^ and the measured velocity broaden- 
ing is 140 km s^^. The derived column density distribution is 
in excellent agreement with the results obtained from model B 
(see Fig. 13. 

In Fig. 13 we compare the ionic column densities predicted 
by model D with the column densities measured by model B. 
The correlation between the measured and predicted column 
density is rather tight, with a measured slope of 1.10 ± 0.07. 
However, one also notes that the measured column densities 
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Fig. 4. Comparison between the ionic column densities calcu- 
lated with the warm model (x-axis is model D) and the ionic 
column densities obtained with the slab model (y-axis is model 
B). The solid line represents the solution where the predicted 
warm and measured slab column densities are equal. 

Table 6. The best fit values for model D, using one warm com- 
ponent and fitting iron separately as a slab component. 



log^ (lO"'' Wm) 



-0.1 
3.5 



(6.3 ± 1.1) X 10^'' m"^ 
(1.6 ±0.3) X lO^** m-^ 



using model B are higher than for model D. This is due to the 
fact that with method B we optimized the outflow velocity. This 
is not possible with model D, where the outflow velocity and 
velocity broadening are tied for all ions. 

3.3. Reliability o1 column density estimates 

In Table [Tsl in appendix A we list all the ions that we fitted 
with the slab model (model B). As the slab model fits both for 
the absorption lines and the corresponding edges, some column 
density determinations could be dominated by an edge mea- 
surement, rather than by absorption lines. Measurements dom- 
inated by an edge, although resulting in quite accurate column 
densities are highly dependent on the correct calibration of the 
instrument and the continuum fit, and thus potentially have a 
large systematic error. A further uncertainty is introduced by 
unknown blending at the longer wavelengths (i.e. above 60 A), 
as the atomic data for lines produced at this wavelength region 
are sometimes poorly known. In Appendix B we discuss those 
ions for which the column density maybe affected by the above 
effects. 

3.4. Narrow emission lines 

Our spectrum shows the presence of a few narrow emission 
lines. Table01ists the measured strength or upper limits of for- 
bidden and intercombination fines of several ions as well as 



Table 7. Narrow emission lines. The Equivalent Width (EW) as 
measured with the LETGS, except for the MgXl, Alxill and 
Si XIII forbidden lines (MEG), and flux are listed. In the last 
column the ionization parameter, ^, where the ion has the high- 
est column density is given. Forbidden lines are indicated by f, 
intercombination lines by i. Below the line the Cv and Ovil 
RRC are listed. The rest wavelengths are taken from Drake 
([1988). 



line 


wavelength 


EW 


flux 


logr 




(A) 


(mA) 


(phm-2 s-i) 


Cvf 


41.435 


283 ± 170 


1.0 ±0.6 


0.1 


Nvif 


29.518 


37 ± 19 


0.26 ±0.13 


0.6 


Nvii 


29.082 


< 17 


<0.1 


0.6 


Oviif 


22.093 


103 ± 14 


0.88 ±0.12 


0.95 


O Vlli 


21.804 


< 5 


< 0.04 


0.95 


Ne ixf 


13.690 


28 ± 10 


0.14 ±0.05 


1.5 


Ne ixi 


13.566 


12 ± 9 


0.05 ± 0.03 


1.5 


Mg Xlf 


9.314 


<4 


<0.02 


1.9 


AlXIlf 


7.864 


4±2 


0.04 ± 0.02 


2.05 


SiXIIlf 


6.739 


5 ± 2.7 


0.03 ± 0.02 


2.2 


Cv 


31.63 


23 ± 16 


0.08 ± 0.06 


0.1 


OviI 


16.77 


51 ±32 


0.08 ± 0.05 


0.95 



Given in 10"^ W m. 



the Radiative Recombination Continua (RRC) observed in the 
spectra. For all narrow emission lines, with the exception of 
the O VII and Ne IX forbidden lines, the wavelength was frozen 
to the value in the restframe of NGC 5548. For the O VII and 
Ne IX forbidden lines we determined blueshifts of —0.009 ± 
0.004 A and -0.005 ± 0.007 A, corresponding to an outflow 
velocity of — 150 km s^^ and —175 km s^^, respectively. These 
blueshifts, however, become negligible if instead, the redshift 
of 0.01717 based upon the 21 cm line (Crenshaw & Kraemer 
[1999) is assumed. 

The intercombination lines as well as the RRC's are unde- 
tectable or very weak. In particular the O VII intercombination 
line is not detected, likely due to blending by the absorption 
fine of Ovi at 21.79 A rest wavelength. The forbidden Ovil 
line seems double-peaked in the MEG spectrum (see online ap- 
pendix C Fig. ll2t . and is poorly fit by a single Gaussian line. 
However, the LETGS spectrum and the earlier MEG spectrum 
of 2000 do not show this line profile. None of the other for- 
bidden lines show any broadening or a double-peak profile. We 
thus conclude that the MEG O VII forbidden line shape is due 
to noise. 



3.5. Fe Ka 

The narrow Fe Ka emission line is clearly seen in the HEG 
spectrum (see online appendix C Fig. 1121) . The flux of the line 
is (0.24 ± 0.08) ph m^^ s^^. Our results are in good agree- 
ment with the EPIC results presented by Pounds et al. J2003> 
and the earlier HEG data presented by Yaqoob et al. (200 D 
(Table|8}. The Fe Ka emission line is also discussed by Yaqoob 
& Padmanabhan ( 



J2004t . There is no evidence for a broadened 
Fe Ka emission line, consistent with the results obtained by 
Yaqoob et al. (.2001,) and Pounds et al. (.2003/ 
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Table 8. Parameters of the Fe Ka line from the present HEG 
spectrum (2002), previous HEG spectrum (2000, Yaqoob et 
al.|2001 1, and EPIC spectrum (2001, Pounds et al. 15005 ). The 
FWHM is given in km s^^, the flux in ph m^^ s~^. 



HEG (2000) 
E (keV) 6.402 ± 0.026 
EW(eV) 133 ±50 
FWHM 4515 ±2650 
flux 0.36 ±0.16 



EPIC (2001) 
6.39 ± 0.02 
60 ± 15 

6500 ± 2200 
0.38 ±0.10 



HEG (2002) 
6.391 ±0.014 

47 ± 15 
1700 ± 1500 
0.24 ± 0.08 




22 22.2 
Wavelength (A) 

Fig. 5. Detail of the LETGS spectrum showing the fit with 
(thick line) and without (thin line) a broad emission line for 
the O VII resonance Une. The profile of the broadened emission 
line is also plotted. 

3.6. Broad emission lines 



Kaastra et al. ( I2002a> detected a 3a significant relativistically 
broadened emission line for N VII Lya and O VIII Lya, and a 
non-re lativistically broadened emission line for C VI Lya. Ogle 
et al. (2D04) detect a broadened C VI Lya line with a FWHM 
of 1100 km s"^ in NGC 4051. No relativistically or non- 
relativistically broadened emission lines were detected in the 
high flux state RGS spectrum (Steenbrugge et al. 20031. The 
LETGS and MEG spectra show clear, broad excess emission 
centered on some of the deeper absorption lines, although none 
show the asymmetric line profile of a relativistically broadened 
line (Fig.|5j. We fit these broad excesses with Gaussians. Due 
to the noise in the data, the broadness of these features, and 
the rather low contrast with the continuum, we fixed the width 
of these Gaussians to 8000 km s^^ (Arav et al. I2002> as ob- 
served for the broadest component of the C IV and Lya broad 
emission lines in the UV. The wavelength was fixed to the rest 
wavelength, i.e. assuming no outflow velocity. Table|9]lists the 
broad emission line fluxes for the LETGS and MEG spectra. 
A detail of the LETGS spectrum of the O VII triplet broadened 
line is shown in Fig.|5] 

To ascertain the significance of these broad emission lines, 
we checked the lines at shorter wavelengths separately in the 
LETGS and MEG spectra. The line fluxes are consistent within 



Table 9. Flux of the broad emission lines, from the simultane- 
ous LETGS, MEG and HEG fit. The wavelength was frozen to 
the rest wavelength of the line, while the FWHM was frozen 
to 8000 km s^^. The lines for which upper limits are detected 
were not used in the further analysis of the data. 



ion 


A 


flux 


EW 




(A) 


(phm-^ s"^) 


(mA) 


Cvf" 


41.421 


1.0 ±0.6 


300 ± 180 


Cv Is2-ls2p ^Pi 


40.268 


<0.6 


< 230 


C VI ls-2p (Lya) 


33.736 


0.5 ± 0.2 


150 ± 80 


N VII ls-2p (Lya) 


24.781 


<0.05 


< 70 


O VII triplet'' 


21.602 


0.56 ±0.13 


130 ± 40 


O VIII ls-2p (Lya) 


18.969 


0.4 ± 0.2 


60 ±30 


Ovil ls^-ls3p ^Pi 


18.627 


0.19 ±0.07 


70 ±30 



" This coincides with the instrumental C-edge, therefore its detection 
is less certain due to possible calibration uncertainties. 

The resonance intercombination and forbidden lines significantly 
overlap, so separate measurements are difficult. 



the error bars, between the MEG and LETGS observations. 
Inclusion of these broad emission lines in the model alters the 
derived column densities for the different ions in the warm ab- 
sorber. Therefore the differences between the model with and 
without a broad emission line is not equal to the flux of this 
line (see Fig. |5Jl. As a result the column densities quoted in 
this paper, which are based on the fit with the broad emission 
lines included, should not be compared directly with column 
densities derived from earlier observations as small continuum 
mismatches are Ukely to remain in those spectra (e.g. Kaastra 
et al. 2002a" Steenbrugge et al. I2003I I in a search for variabil- 
ity of the warm absorber. As an example, the logarithm of the 
derived column density for O VII (using one slab component) 
in a model without a broad emission line is 21.72 m~^, com- 
pared to 22.18 m^^ if the broad emission lines are included. 
The broad emission lines are more easily detected than in pre- 
vious observations, this partly results from the low continuum 
flux level, and thus the larger contrast during the present obser- 
vation. Further discussion on these broad emission lines will be 
given in a forthcoming paper (Steenbrugge et al. l2004> . 



4. Timing analysis 

Kaastra et al. (2004) give a detailed analysis of the lightcurve 
in different energy bands. Here we focus on the possible differ- 
ences in the warm absorber due to flux variations. 

4. 1. Sliort term variability 

First we studied possible differences in the HETG and LETG 
spectrum, as the average flux level differed by 30 % between 
both observations on a timescale of about 300 ks. In our first 
approach, we have taken the fluxed MEG spectrum and folded 
this through the LETGS response matrix. This is possible be- 
cause of the two times higher spectral resolution of the MEG as 
compared to the LETGS. We adjusted the MEG continuum to 
match the higher LETGS continuum. Restricting our analysis 
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Fig. 6. The transmission for the best fit model with three xabs 
components for the LETGS (thin sohd Hne) and the HETG 
(dotted Hne) spectra. Note that the difference is largest, 10 % 
in continuum transmission and 35 % for the largest difference 
in absorption lines, for lowly ionized ions which have absorp- 
tion lines around 20 A. Here the MEG effective area is already 
rather low, and the errors on this lowest ionization parameter 
are large. 

to the 1 .5-24 A band, where both instruments have sufficient 
sensitivity, we find good agreement between both spectra. 

We find no significant difference in the equivalent width of 
the strongest absorption lines such as the Lya and Ly f3 lines 
of Ne X and O VIII and the resonance lines of Ne IX and O VII. 
There may be a small enhancement of the most highly ionized 
lines in our spectrum (Si XIV and Mg XII Lya and the Si XIII 
resonance line), corresponding to a ^50 % increase in ionic 
column density in the LETGS observation. For the Si XIII line 
this is uncertain because the line coincides with an instrumen- 
tal edge in the MEG spectrum. Moreover, at the high ionization 
parameter where these ions are formed (log^ ~ 2.5) the opac- 
ity of the O VIII Lya line is even higher than for the Mg and 
Si lines. Since we do not observe an enhanced column density 
for O VIII, we conclude that the evidence for enhanced opacity 
at high ionization parameter during the LETGS observation is 
uncertain. A formal fit to the difference spectrum shows that 
any additional high ionization component emerging during the 
LETGS observation has a hydrogen column density less than 
10^'' m^^, or an order of magnitude less than the persistent 
outflow. 

In a different approach, we have fitted the warm absorber in 
the LETGS spectrum between 1 .5 and 24 Ausing a combina- 
tion of three xabs components, the results are listed in TableFTol 
Fitting the HETG spectra with the same ionization parameters 
and column densities, but fixing the continuum parameters to 
the best fit parameters (see Table |2jl we find an excellent fit. 
Allowing the ionization and column densities to be free param- 
eters, we obtain the results listed in TablellOl The improvement 
for allowing the ionization and column densities to be free is 
Ax^ = 10 for 4242 degrees of freedom. The maximum differ- 
ence between the two data sets occurs for the lowest ionization 



Table 10. Comparison of the best fit parameters for the LETGS 
and HETGS spectra between 1 .5 and 24 A. The labeling of the 
different components follows Kaastra et al. J2002a> . 





LETGS 


HETGS 


LETGS 


HETGS 




logr 


logr 


log Nh' 


log Nh' 


1 


2.26 ± 0.09 


2.3 ± 0.3 


25.47 ± 0.06 


25.36 ±0.08 


2 


1.77 ± 0.04 


1.9 ±0.1 


24.85 ± 0.14 


25.02 ±0.12 


3 


-0.2 ± 0.2 


-0.1 ±2.5 


24.33 ±0.16 


24.8 ± 0.4 



" Ionization parameter in 10 W m. 
Hydrogen column density in m"^. 



component 3. For the HETG spectra this component has log 
f = —0.1 ± 2.5, this clearly indicates the lack of sensitivity 
at the longer wavelengths in the HETG spectra. Fig. |6] shows 
the transmission for the best fit model for LETGS (solid line) 
and the HETGS (dotted line). Note that the largest discrepancy 
is around 20 A, where the effective area of the MEG is rather 
low, and the lines are for the lowest ionized component 3. The 
flux in the soft band increased by ^50 % between the HETGS 
and LETGS observation. Therefore, an increase of ^ by 50 % is 
excluded for component 2. However, our data cannot rule out 
that the ionization parameter responds linearly to continuum 
flux enhancements for components 1 and 3. We conclude that 
the difference in ionization parameter of the warm absorber be- 
tween the HETG and the LETGS observation is small and be- 
low the sensitivity in log ^ of ~0.15. 

A similar situation holds for the time variability during the 
LETGS observations. Despite the significant flux increase dur- 
ing the observation (see Kaastra et al. 2004 1, the signal to noise 
ratio of our data is insufficient to rule out or confirm a signif- 
icant response of the warm absorber to the change in ionizing 
flux. 

4.2. Long term spectral variations 

We have tested for the presence of long term variations in the 
warm absorber of NGC 5548 by comparing our best fit model 
for the present data with the LETGS spectrum of December 
1999 (Kaastra et al. 2000 2002a i. To do this, we have fit the 
1999 spectrum keeping all parameters of the warm absorber 
and all fluxes of narrow or broadened emission lines fixed to 
the values of our 2002 spectrum. Only the parameters of the 
power law and modified blackbody component were allowed to 
vary. Then we subtracted the relative fit residuals for our 2002 
spectrum from the fit residuals of the 1999 spectrum (Fig.0. 
The relative fit residuals r are defined by r — {o — ra)/m 
where o is the observed spectrum, m the best fit model spec- 
trum. Subtracting the fit residuals instead of dividing the spec- 
tra allows us to take out continuum variations. Further, as the 
higher orders are not subtracted any changes in the continuum 
will produce broad band features in the ratio of two spectra. In 
Fig.0we normalized the difference of the fit residuals by the 
standard deviation of this quantity. This erases all remaining er- 
rors due to either shortcomings in the effective area calibration 
or in the spectral modeling. 

The figure shows clear evidence for five features with more 
than 2.5(7 significance. There is evidence for an increasing flux 
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Fig. 7. Fit residuals of the 1999 LETGS spectrum minus the fit 
residuals of the 2002 spectrum, normalized as described in the 
text. For clarity, the residuals have been rebinned by a factor of 
8. The dashed lines indicate the 2.5 tr significance level used in 
the analysis. 

Table 11. Features in the difference spectrum of the 1999 and 
2002 LETGS data. We list Gaussian centroids A (in the rest- 
frame of NGC 5548), cr and peak value p; peaks are expressed 
as a fraction of the 2002 spectrum at the same energy. 



A (A) 


a (A) 


P 


16.10±0.06 


10+""'-' 

"■^"-0.04 


-0.25±0.09 


22.36±0.02 


0.03±0.02 


+0.85±0.29 


27.18±0.07 




-0.27±0.12 


34.13±0.08 


0.12±0.07 


+0.34±0.14 



toward smaller wavelengths starting around 5 A and resulting 
in about a 20 % excess at 2 A in the 1999 spectrum. This cor- 
responds to a difference in the continuum and not in the warm 
absorber. The remaining features are not due to continuum vari- 
ations (Table [m . As the optical and UV broad emission lines 
are known to be variable on timescales as short as a day, we 
identify the broader features with changes in the broad emis- 
sion lines. We identify the 34.13 A feature as the red wing of 
a broad C VI Lya line (restframe wavelength 33.74 A, hence 
velocity of +3500 ± 700 km s'^). Similarly, the 27.18 A fea- 
ture can be identified as the red wing of C VI Ly7 (restframe 
wavelength 26.99 A, hence velocity of +2100 ± 800 km s^^). 
Apparently, the red wing of the C VI Lya line has decreased 
while the red wing of Ly7 increased and no change in the Ly/3 
line. However, the process that would cause the Lya line to 
decrease while the Ly7 line increases is unknown. The intrin- 
sic width, cr, of the variable parts of both lines corresponds to 
about 900-1000 km s^^, this is much smaller than the assumed 
FWHMof8000kms-i. The 16.10 A feature can be identified 
by a red wing of Ovill Ly/3 (restframe wavelength 16.01 A, 
hence velocity of +1700 ± 1100 km s^^). A change in the 
strength of the broad emission line due to the O VII triplet can 
also explain the dip at 21 — 22 A. 



Finally, the 22.36 A feature most likely has a different ori- 
gin. First, it is intrinsically naiTow and poorly resolved (the up- 
per limit to its Gaussian width corresponds to 700 km s^^). 
Its wavelength coincides within the eiTor bars with the wave- 
length of the strongest O V absorption Hne (22.374+0.003 A). 
The 2002 spectrum shows a deep O V absorption line, which 
is much weaker in the 1999 spectrum. This demonstrates that 
the Ov column is variable on the time scale of a few years. 
Note that the O V line is by far the strongest absorption line 
in the LETGS spectrum for the ionization parameter log ^ = 
—0.2, the ionization parameter where O V reaches it maximum 
concentration. The Fe M-shell lines are formed at the same ion- 
ization parameter, but are much weaker 

5. Discussion 

5. 1. Comparison between the warm absorber 
models 

In our model B we determined the column density, the aver- 
age outflow velocity, and the velocity broadening for all ob- 
served ions separately (Table|3i. The total column densities de- 
termined using this method, for the six ions C VI, O V, O VII, 
Ovill, Neix, and Sixill are consistent with the results ob- 
tained using method A, with fixed velocity structure. This gives 
confidence that the column densities as listed in Tables |4] and 
I15lare reliable, even if there is substructure to the lines. The av- 
erage outflow velocity determined for these six ions is consis- 
tent with the outflow velocity of the deepest UV components. 
Comparing our measured column densities with those obtained 
in the UV for the five outflow velocity components, we con- 
clude that the column density measured for the —160 km s^^ 
component is in fact the blend of UV components 4 and 5. Our 
—530 km s^^ component is a blend of UV components 2 and 
3. Components 2 through 5 form one unresolvable blend in X- 
rays, and only component 1 is clearly separated. 

5.2. Outflow velocity 

Our method B does not resolve the full velocity structure of 
the outflow. Therefore one should take care in interpreting the 
relation between the ionization parameter and the average out- 
flow velocity (Fig. O. Namely, all line profiles probably are 
a blend of five outflow velocities, but due to limited signal to 
noise ratio and spectral resolution, components 3 and 4 of the 
UV dominate the blend. 

In general, there is a trend that more highly ionized ions 
have a larger average outflow velocity. This could indicate that 
the —1040 km s^^ outflow velocity component becomes more 
prominent for the more highly ionized ions. The —530 km s^^ 
outflow velocity component dominates for lower ionized ions 
observed in X-rays and the UV band. To further study this ef- 
fect we plot the total column density as a function of outflow 
velocity (Fig.|8j for the six ions of Table|3l We also added C IV 
and N V taken from Crenshaw et al. ( 2003 1 and O VI taken from 
the non-simultaneous FUSE data by Brotherton et al. J2002> . 
For the UV data we added the two components with the low- 
est outflow velocity to compare it with the —160 km s^^ X- 
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Fig. 8. The ionic column density versus the logarithm of the 
ionization parameter for the three outflow velocity components 
as measured for six ions in X-rays, 2 ions (C IV and N v) mea- 
sured simultaneously with HST STIS (Crenshaw et al. I2003> 
and the lower limits for Ovi from non-simultaneous FUSE 
data (Brotherton et al. 120021 taking their preferred uncovered 
model; Arav et al. I2003> . In calculating the column densi- 
ties for the UV data we kept the highest outflow value sepa- 
rate, added the two middle ones and the two lower ones (see 
Sect. 15. 2> . For O VII the ionization parameter was decreased by 
0.05 for easier identification. The open circle indicates the total 
column density for O VI obtained from the LETGS and MEG 
spectra. 



ray component. The two middle outflow velocity components 
were added to represent the —530 km s^^ component. Finally, 
for O VI we indicate the total column density as measured in 
the X-rays. From Fig. |8]we note that the ionic column den- 
sity of the —1040 km s^^ velocity component is the smallest 
of the five velocity components for low ^, while it becomes the 
largest for highly ionized gas. There is thus a clear difference in 
ionization structure between the —1040 km s^^ outflow com- 
ponent and the lower outflow velocity components. A possible 
explanation for this difference in ionization structure is that the 
outflowing wind at —1040 km s^^ is less dense, leading to a 
higher overall ionization of the gas, while overall the column 
density is smallest. 

5.3. Velocity broadening 

In Table|3]we list the measured Gaussian velocity broadening a 
of the three velocity components as well as the velocity broad- 
ening of the UV components 1,2 — 4, and 5. For components 
2 — 4 we find a slightly smaller X-ray width, while for com- 
ponent 5 we find a larger width as compared to the UV lines. 
We attribute this to the blending of (a part of) component 4 into 
component 5 . We note that U V component 4 has ct = 62 km s ^ ^ 
compared to 68 km s^^ in component 3 and only 18 km s^^ in 
component 2, hence we need approximately half of component 
4 to blend into component 5 in order to explain the difference. 



Interestingly, we find a much smaller velocity broadening 
for UV component 1, the only one we can resolve. Modeling 
the O VIII Lya line with a column density of 10^^^ m~^ and a 
velocity dispersion cr of 40 km s~ ^ , we find that the line is heav- 
ily saturated, and effectively produces a FWHM of 240 km s^ ^ . 
This is similar to the 222 ±18 km s~^ value measured from 
the UV spectrum. We are, however, able to obtain the velocity 
broadening a from the line ratios of the non-saturated Lyman 
series of O VIII and the other five ions. 

In Table |4] the velocity broadening a is listed based upon 
model B. It is small compared to the velocity range of the 
components measured from the UV spectra. This indicates that 
we do not fully resolve the blend of the five outflow veloc- 
ity components. The velocity broadening in our case is mainly 
determined from equivalent width ratios of absorption lines 
from the same ion. The velocity broadening for the majority 
of absorption lines is less than 220 km s^^. This is consistent 
with detecting mainly the —667 km s^^, —530 km s^^ and 
—336 km s^^ outflow components, which are the dominant ve- 
locity components in the UV. These three velocity components 
form an unresolved blend in the X-rays. Plotting the velocity 
broadening versus the ionization parameter produces a scatter 
plot. 

6. Ionization structure 

An important question about the physical nature of these out- 
flows is whether they occur as clumps in an outflow or in a 
more homogeneous wind. In the first case these clumps should 
be in pressure equilibrium with the outflowing less dense wind, 
with a finite number of solutions for temperature and ioniza- 
tion ^ for the same constant pressure ionization S value. In the 
second case one expects a continuous ionization distribution, 
which depends on the density of the particular part of the wind. 
Thus less dense material will be higher ionized, while more 
dense regions will be lower ionized. In this section we elaborate 
on both scenario's, confronting them with our observations. 

6. 1. Clumps in a wind 

In several AGN outflow scenario's, instabilities in the wind 
may cause the formation of clumps. Such clumps can survive 
for longer times if they are in pressure equilibrium with the 
suiTounding wind and are on the stable branch of the T(S) 
curve. Where the ionization parameter for constant pressure, 
S, is given by S = L/{Aircr'^P) = 0.961 x lO'^ ^/T, with L 
the luminosity, P the pressure, r the distance from the ionizing 
source, c the speed of light and T the temperature. Ionization 
codes, like XSTAR (Kallman & Krolik 1999) predict a spe- 
cific ionization range for which pressure equiUbrium can occur 
(Krolik & Ki-isslODD. Ki-ongold et al. (20031 claim that there 
are only two ionization components in such an equilibrium in 
the Seyfert 1 galaxy NGC 3783. Netzer et al. (2003) needed 
three ionization components in pressure equilibrium to fit the 
same NGC 3783 spectrum. Ogle et al. ( 2004 i find that in the 
case of NGC 4051 the different ionization components are not 
in pressure equilibrium. Another possibility is that the clumps 
are magnetically confined. 
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Table 12. Fits with a finite number of ionization components 
to the column densities of all elements but iron derived with 
model B, assuming solar abundances. The number of ionization 
components is indicated by N. In the last column we give the 
significance according to an F-test of the added component. 



N x'/d.o.f. log C 



10"'' Wm 



logTVn 



sign. 



1 147/22 

2 86/20 

3 83/18 



1.2 

2.45, 

2.47, 



1.15 
1.17, 



-0.70 



25.5 

26.1, 25.4 
26.1, 25.4, 23.( 



89 

53 



Table 13. Fits with a finite number of ionization components 
to the iron column densities for the RGS data (Steenbrugge et 
al. l2003> . The number of ionization components is indicated by 
N. In the last column we give the significance according to an 
F-test of the added component. 



N 


X^/d.o.f. 


log? 


log A^H 


sign. 






10-^ Wm 


m-2 




1 


79/15 


0.22 


25.2 




2 


22.9/13 


2.50, 0.20 


25.7, 25.2 


98.8 


3 


10.5/11 


2.57, 1.61, 0.20 


25.7, 25.1, 25.2 


90.9 


4 


7.2/9 


2.57, 1.62, 0.65 


25.7, 25.1, 25.0 


72.7 






-0.03 


25.0 




5 


4.9/7 


3.10, 2.48, 1.56 


25.8, 25.6, 25.1 


71.1 






0.65, -0.03 


25.0, 25.0 





In order to test the presence of a finite number of ionization 
components, each with its own value for ^ and A^h. we fitted 
the ionic column densities obtained with model B (Table FTSt to 
a model with a finite number of ionization components. Since 
there are potential problems with the iron ionization balance 
(see Sect. 13 .2. 3t . we did our analysis separately for iron and the 
other elements. The results of these fits are listed in Table 1 121 
and Tableini 

Using the measured column densities for all ions, except 
iron, we notice that the fit does not improve if more than three 
ionization components are fitted. The program even prefers to 
split up one component rather than add another in the case we 
fitted for five ionization components. The fit is never statisti- 
cally acceptable, which is unlikely due to abundance effects as 
most elements have abundances consistent with solar. As an ex- 
tra test we decided to fit oxygen, silicon and sulfur separately. 
However, due to the smaller span in ionization range and the 
fewer points (maximum of five), all were well fit with three 
ionization components. 

For iron we decided to use the column densities measured 
by the RGS (Steenbrugge et al. I2003> . as here several more 
ions have a measured column density and the eiTors are some- 
what smaller For Fexil and Fexxiv, which have the same 
measured column density in the RGS and Chandra spectra, we 
used the smaller Chandra error bars. To fit these column den- 
sities adequately, one needs at least five ionization components 
(log ^ = 3.10, 2.48, 1 .56, 0.65, -0.03). In Fig.|9lwe show the fit 
to the measured column densities with four and five ionization 
components. The main problem with the fit for four ionization 
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Fig. 9. The column densities for iron as measured by the RGS 
instrument, with the fit for a model with four (dash-dot line) 
and five (thick solid line) ionization components. The fit with 
four ionization components clearly underpredict the Fe XXIV 
column density. 



components is Fe XXIV, which is severely under-predicted. As 
the RGS and Chandra spectra measure the same column den- 
sity for this ion, this is a rather certain measurement. 

We conclude that we need at least five different ionization 
components to explain the measured ionic column densities in 
NGC 5548. These components span a wide range in ionization 
parameter (log ^ between —0.03 and 3.10). Are these ioniza- 
tion components in pressure equilibrium? In Fig. [TO] the ion- 
ization parameter if pressure is held constant (S) versus tem- 
perature is shown for the SED used in the current analysis and 
the earlier RGS and LETGS analysis. The marginally stable 
part of the curve indicates a constant nT value, therefore com- 
ponents with a different ^ are in pressure equilibrium on the 
marginally stable branch from log ^ = 1.3 to 2.7. The curve 
has nowhere a negative slope, the smallest positive slope is be- 
tween log — 1.4 — 1.5. It should be noted that the difference 
between the results produced by different codes can be large, 
but substituting CLOUDY 95.6 (Ferland'2002^ for XSTAR we 
find that still not all components are in pressure equilibrium. 

For the three different ionization components as observed 
in the earlier LETGS and RGS spectra of NGC 5548 (Kaastra et 
al. ,2002a. Steenbrugge et al. .2003j . the lowest ionization com- 
ponent has log ^ = 0.4, which cannot be in equilibrium with 
the other two components at log ^ = 1.98 and 2.67. If we take 
the five ionization components from the iron analysis, then the 
highest and the lowest ionization component cannot be in pres- 
sure equilibrium. This clearly suggests that the current photo- 
ionization balance models are either too simple, or that at least 
the lowest ionization component is confined by another pro- 
cess, possibly magnetic. As the lowest ionization component 
does have the same kinematics as the higher ionized absorbers, 
all absorbers must co-exist in a single confined outflow. The 
result that not all ionization components are in pressure equi- 
librium is rather robust, even accounting for inaccuracies in the 
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Fig. 10. The temperature versus ionization parameter for con- 
stant pressure S diagram for the SED used in the present analy- 
sis (thick hne) and for the SED used in earlier LETGS (Kaastra 
et al. l2002a> and RGS (Steenbrugge et al. l2003l analysis (thin 
line). The values corresponding to the currently used SED 
are indicated on the x-axis on top. The five components derived 
from fitting iron in the RGS spectrum are indicated. The dotted 
lines indicate the boundaries of the marginally stable branch 
for the SED used in analyzing the RGS data; the dash-dotted 
lines are the same boundaries, but for the current used SED. 



photoionization balance codes, in particular for lowly ionized 
iron, or the inclusion of the blue bump. 

Another possible method to distinguish between a contin- 
uous and discrete ionization is to look for spectral variability. 
No spectral variability of the warm absorber, with the excep- 
tion of O V, was detected although the source was observed at a 
2 - 1 keV lu minosity of 4.9 x 10^6 W (1999 LETGS, Kaastra 
et al. I2002a> . 5.7 x lO^^ W (2001 RGS, Steenbrugge et al. 
l2003t and 4.0 x 10^^ W (present LETGS spectrum). However, 
if the ionization parameter varies linearly with luminosity, then 
no variation in ionization is expected to be detectable. Further 
complicating the possible detection of changes in the ioniza- 
tion parameters is the fact that in the UV the different velocity 
components show opposite column density variations for the 
same ion. The net effect will be hard to discern in the X-rays 
for velocity components 2 — 5. Component 1, which can be 
resolved with Chandra is the most variable component, and 
should be studied in the future for signs of variability. A fur- 
ther lack of variation in the ionization structure of component 1 
possibly indicates that the ionization distribution is continuous. 
In the case of two or three ionization components, a luminosity 
change should result in a shift in the ionization parameters. In a 
continuous model one would expect a decrease in column den- 
sity for the lowest ionization states and an increase in highly 
ionized material for an increase in luminosity. However, as the 
warm absorber is detected up to Fe XXIV, and our sensitivity for 
Fe XXV and Fe XXVI is insufficient, this is hard to measure with 
this data set. However, for the long NGC 3783 XMM-Newton 
observation this behaviour was indeed detected. Behar et al. 
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Fig. 11. The total hydrogen column density A^h, assuming so- 
lar abundances (Anders & Grevesse ll989> derived using Eq.|3l 
plotted versus ionization parameter. The ionic column densi- 
ties were taken from Tablel^and Tablell5lassuming a velocity 
broadening of 140 km s^^. For clarity, no upper limits have 
been plotted. The best fit results for model D are plotted as the 
two crosses connected by a dotted line. 



(I2003t studying the RGS spectra detected no spectral variabil- 
ity in the lowly ionized absorber, and in particular no change 
in the iron UTA. Reeves et al. (2004) however did detect spec- 
tral variability of the highly ionized absorber with the EPIC pn 
instrument. 

6.2. Continuous ionization distribution 

We decided to derive the total hydrogen column density as a 



function of ^ for each ion, similar to Steenbrugge et al. J2003> . 
This in contrast to the warm model where we determined the 
total column density at the two extrema. Contrary to the pre- 
vious analysis we did not assume that the ion column densi- 
ties are predominantly determined by the ionization parameter 
at which the relative concentration is at its maximum. In the 
present analysis we took into account that each ion is formed 
over a range of ionization parameter similar to the warm model. 
We implemented this by assuming a power-law distribution, as 
this function well describes the RGS data (Steenbrugge et al. 
2003 1, and appears to be sufficient in the warm model: 



(1) 



with a given normalization A and index a. We calculate 
the predicted ionic column densities A'i p for this distribution 
similar to the warm model. As the ionization parameter range 
in reality will be finite, this power law needs to have a cut-off 
value, which was chosen to be log ^ = 4. We will use the results 
from model B to determine A and a. If .fi(0 is the ionic con- 
centration relative to hydrogen at a given ionization parameter 
^, then Nup is given by N^p = fiiOM"'^^ d^. For each 
ion the effective ionization parameter is defined as 
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(2) 
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This value of is typically 25 % larger than the value of 
^ for which reaches its maximum. This increase is solely 
due to our assumption of a power-law distribution. This can be 
used to derive the equivalent hydrogen column density Nn ifor 
each ion (not to be confused with the integrated column den- 
sity over all ^ values) from the observed ionic column densi- 
ties M obs, since obviously we should have iVj.obs = M.p- This 
yields equation|3j 



diVH(6: 



i.obs 



(3) 



The calculated total hydrogen column density using this 
method, instead of assuming that the ion is formed at the ion- 
ization parameter where /i(^) has its maximum, is larger by at 
most a factor of 3. Ion concentrations were determined using 
a set of runs with XSTAR (Kallman & Krolik 1999l for a thin 
photoionized layer with solar abundances (Anders & Grevesse 
[T9891 . 

The derived hydrogen column densities A^h.i versus the 
ionization parameter are presented in Fig. ^2 For lower ion- 
ized ions (log ^ < —0.7) we can only determine upper limits 
to the column density, in contrast to the earlier RGS results 



(Steenbrugge et al. I2003t . This results from the lower effec- 
tive area of the LETGS compared to the RGS in the iron UTA 
region. However, due to the longer wavelength band, we do de- 
tect many previously unobservable medium ionized ions. This 
strengthens the relationship between the hydrogen column den- 
sity and ionization parameter in the range log ^ = to 1 . The 
column density increases by nearly 2 decades for a 3 decades 



increase in ^, consistent with Steenbrugge et al. J2003> . 

A formal fit of the column densities of all ions except iron 
yields A = 1024-4±0.i (^t log ^ = 0) and a = 0.40 ± 0.05 with 

= 107 for 30 d.o.f. A similar fit to the iron ions only yields 
A = l025 0±0-2 (at log C = 0) and a = 0.2 ± 0.1, with a = 16 
for 13 d.o.f. 

The relationship obtained is rather tight, and the assump- 
tion of solar abundances seems justified, if we do not include 
those ions which may suffer from systematic bias in their col- 
umn density determinations ("Sect l3.3l l. If the same analysis is 
repeated for sodium, aluminum, argon and nickel, we find that 
these elements are overabundant compared to solar, except for 
sodium. This strengthens our conclusion that these ions may 
suffer from strong systematic errors in their column density es- 
timates. These ions are therefore not plotted in Fig. II II 

7. The outflow 

7.1. Outflow geometry 

In trying to get a physical picture of the outflows in AGN it is 
important that we know the geometry, as well as the ionization 
structure. By estimating the opening angle of the outflow, we 
can distinguish between a spherical outflow model or a model 
with localized streams. The model by Elvis ("20001 prefers a 
localized stream. The model by Urry & Padovani (.1995j . on 
the other hand predicts a less coUimated outflow. 

Similar to Steenbrugge et al. f2003) we can obtain an up- 
per limit to the opening angle for the outflow (eq.|3. We use 



Table 14. The upper limit to the opening angle as calculated 
from eq.l^for the five outflow velocities measured in the UV, 
over the ionization range observed in the Chandra spectra. The 
angles are given in sr. 
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4 
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2 




1 


log r 


-166 


-336 




-530 


-667 




-1041 




km s^^ 


km s^^ 




km s^^ 


km s^ 


1 


km s""^ 





7x10-"* 


3.5x10^ 


'4, 


2x10"" 


2x10" 


-4 


1x10"* 


1 


7x10"^ 


3.5x10" 


-3 
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2x10' 
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1x10"^ 
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3.5x10" 


-2 
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1x10"^ 
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7x10"^ 


3.5x10" 


-1 


2x10"^ 


2x10 


-1 


1x10"^ 



" in 10"'-* W m. 

the mass conservation formula (Mjoss = ^iTipTierQv), with tq 
the distance from the ionizing source. He the electron density, 
nip the proton mass, v the outflow velocity and fl the open- 
ing angle. Further we assume that the system is stationary and 
that the mass loss rate is equal or less than the mass accretion 
rate Mace = L/rjc^, with L the bolometric luminosity and 
the accretion efficiency. The outflow velocities of the wind are 
taken from the UV measurement. This leads to the following 
constraint on the solid angle Vl: 



n < 

where we used: 

L 



Tier" 



(4) 



(5) 



We assume a Schwarzschild black hole, with an accretion 
efficiency 77 = 0.057 (Thorne 1974'l. For a Kerr black hole, with 
an efficiency of 0.30 (Thorne 1974) , the opening angle is re- 
duced by a factor of 5.4. Even higher efficiencies will further 
reduce the opening angle. Assuming that only 50 % of the mass 
that is accreted is lost through an outflow halves the upper limit 
to the opening angle. 

As the ionization parameter, as well as the outflow veloc- 
ities measured in this analysis closely resemble those found 
with the RGS, also the resulting upper limits for 51 have a sim- 
ilar range (but now corrected for the extra 4tt). Table \T4\ lists 
the upper limits to the opening angle for the five different out- 
flow velocity components observed in the UV and the ioniza- 
tion range observed in the X-rays. 

We thus find that the opening angle of the outflow is very 
small, down to lO^** sr for the lowest ionization component. 
This implies that the wind occurs in the form of narrow stream- 
ers. The range in ionization parameters at a given outflow ve- 
locity component then implies density stratification perpendic- 
ular to the outflow. 

A further quantitative diagnostic of the thickness of the out- 
flow can be derived from Fig.[^ There is a clear power-law re- 
lation between the total hydrogen column density and the ion- 
ization parameter This power-law relation is also established 
by the fact that one can fit the spectra with the warm model, 
using only two points, the lowest and highest ionization state 
observed. From Fig. we find a power-law slope a = 0.40 
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± 0.05, and a normalization log A = 24.4 ±0.1 at log ^ = 
(see Sect. 16. 2> . Similar to Steenbrugge et al. |2003 1 we com- 
bine Eq.^with Eq.|5]and use dA^H = dr with r the distance 
from the axis of the streamer, resulting in a simple differential 
equation for d^/dr with formal solution n/rio = (1 + r/vo)^'^ 
where Uo is the density at the axis and Tq = A£^"^^/no{a + 1) 
and /3 = l/{a + 1). As we do not know the distance Tq of the ab- 
sorber to the nucleus of the AGN, we do not know Uo', however 
weknowrono= 1.0 x lO^^ m'^ at ^ = and /3 = 0.71 ±0.01. 
Accordingly, at large distance from the axis of the streamer 
(low density, high ionization parameter), the density scales as 
n |,-o.7i Steenbrugge et al. (2003 1 found a similar result, 
but due to the detection and overabundance of the iron UTA 
with the RGS, the power-law slope was less well determined. 

7.2. Mass loss through outflow 

Outflows from AGN are potentially important contributors to 
the enrichment of the IGM. A critical factor is the amount of 
mass that escapes from the host galaxy through the observed 
outflows. Here we estimate an upper limit to the mass loss and 
discuss briefly the difficulties associated with this estimate. 

First the outflowing gas has to escape the potential of the 
massive black hole. Matter can escape if its radial velocity Vr 
is larger than the escape velocity = \/2GM /tq. Inserting a 
black hole mass of 6.8 x 10^ Mq (Wandel 2002) we find that 
matter escapes only if ro > (0.58/u^) pc with Vr expressed in 
units of 1000 km/s. The measured outflow velocities are only 
lower limits, namely the velocity component along our line of 
sight, and range between —1041 kms^^ and —166 kms^^. 
Taking the measured outflow velocities, then tq needed for the 
outflow to escape corresponds to 0.5-20 pc. The location of the 
wind is thus important in determining whether the gas will es- 
cape or not. Crenshaw et al. (2003) argue that the wind should 
cover at least a part of the inner, high ionization narrow line re- 
gion (NLR), while Arav et al. (2002) argue that the narrow line 
region is not covered. The NLR spans at least the distance range 
of 1 pc (highly ionized) to 70 pc (lowly ionized) (Kraemer et 
al. ll998> . If the argument of Crenshaw et al. (2003) holds, then 
in particular the highest velocity gas may be able to escape. If 
the Vr is significantly higher than the measured v, the lower 
velocity components will also escape. 

After escaping the immediate environment of the central 
black hole, the gas must escape from the bulge and halo 
of the galaxy. The stellar velocity dispersion of the bulge is 
180 kms^^ (Ferrarese et al. i200 1 1, and hence it is likely that 
most of the gas can escape the bulge provided it has not been 
decelerated. The effects of the disk and halo of NGC 5548 are 
expected to be of equal or less importance. 

We neglect the possibility that ram pressure in the host 
galaxy terminates the outflow. Wilson & Ulvestad J1982> note 
the possibility that the linear radio structure observed in the 
host galaxy of NGC 5548 is due to material ejected from the 
nucleus and stopped due to ram pressure. This seems to indi- 
cate that the outflows as observed in the UV and the X-rays can 
escape the nucleus, whether or not the host galaxy. 



If all gas escapes, then our earlier argument that the total 
mass loss through the wind must be smaller than the accretion 
rate through the disk leads to an upper limit of the mass loss 
of about 0.3 Mq yr^^. NGC 5548 had a major interaction with 
another galaxy some 0.6-1.0 Gyr ago (Tyson et al. [T998i; as- 
suming that the AGN phase lasts for at least 0.6 Gyr with a 
steady mass loss rate, the upper limit to the mass enrichment of 
the IGM is 2 X 10^ M©. 

8. Summary and conclusions 

Our long Chandra observation of NGC 5548 allowed us to ob- 
tain a unique, high signal to noise, high spectral resolution, 
broadband spectrum of this Seyfert 1 galaxy. 

The fluxes measured for the Fe-K line and the narrow emis- 
sion lines are within 1 a consistent with earlier Chandra and 
XMM-Newton spectra. There is thus no evidence for long-term 
variability. The only narrow feature that shows column density 
variations is O V. 

We found clear evidence for broad emission lines, similar to 
the optical and UV Broad Line Region lines. In particular C VI, 
O VII and O VIII showed significant broad lines. From a com- 
parison of our spectrum taken in 2002 with the earlier LETGS 
observation of 1999, we find significant changes in the red wing 
(at 1700-3500 km s^^) of C VI lines and the O VIII Ly/3 line. 
A more extensive discussion of these broad lines will be given 
in a forthcoming paper (Steenbrugge et al. 2004). The presence 
of these broad lines affects our estimates of column densities 
in the warm wind, by factors up to 3, a similar situation as has 
been found in the UV band (Arav et al. 2002). 

We find evidence for the presence of an inner-shell X-ray 
absorption line of N V at 29.42 A. The long exposure time and 
the combination of LETGS and HETGS spectra allows us to fit 
the warm absorber using three different spectral models. Based 
upon these models, we conclude that there is a good agreement 
between the properties of the outflow as measured through UV 
absorption lines and X-ray absorption lines, although in the X- 
ray band we do not fully resolve the spectral lines. But line 
centroids and derived line width are in good agreement. We 
find that the highest velocity component 1 at —1040 km s^^ 
has a different ionization structure than the other components 
2-5. 

We have compared our results with models for the ioniza- 
tion structure of a wind. We conclude that our data are not in 
agreement with a model with discrete clumps that are in gas 
pressure equilibrium with the surrounding medium. This con- 
clusion is based upon the following arguments: 

1 . the need for at least five discrete ionization components to 
span the broad range (at least 3 orders of magnitude) of 
ionization parameter ^; 

2. we find, assuming our SED, a limited range in ,^ for which 
multiple solutions of constant S can co-exist. The lowest 
ionization parameter determined from the data is incom- 
patible with this range; 

3. the fact that there are five discrete velocity compo- 
nents, instead of a continuous velocity range. At least the 
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— 1040 km s component spans the observed range in ion- 
ization. 

Our data are in agreement with a model with a continuous 
distribution of iVn as a function of Comparing the mass out- 
flow through the wind with the accretion rate onto the black 
hole, we derived upper limits to the solid angle sustained by 
the outflowing wind. These upper limits are of the order of 
10"'' sr for the lowest ionized gas in velocity component 1. We 
conclude that the wind occurs in the form of narrow streams 
with density variations perpendicular to the flow velocity. Most 
likely these streamers are caused by accretion disk instabilities 
and cross our line of sight when being radiatively accelerated. 
Assuming density stratification, the dependency of A^h upon 
gives the density profile across these streamers. At large dis- 
tances the density is proportional to r"°'^^, with r the distance 
from the densest part of the streamer. 

Unfortunately the LETGS is not sensitive enough to de- 
tect significant variations of the warm absorber as a response to 
the continuum flare occurring during the LETGS observation. 
Reverberation studies with more sensitive X-ray instruments 
will allow us to determine the densities and hence the location 
of the wind. However regular monitoring of NGC 5548 in both 
X-rays and UV will allow us to derive important constraints on 
the long term variability of the wind. A clear example is the 
variability of O V between 1999 and 2002 that we found with 
our observations. 
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Table 15. The best fit column densities (in m^^) as measured using model B. The outflow velocity (in km s^^) was taken from 
Table 0] or a function of ionization and frozen. The velocity broadening (in km s^^) was also frozen during the fit. In the last 
two columns we list the ionization parameter (in 10^^ W m) for which the ion has its maximum column density. The column 
densities quoted are for those measured in the fit including the broad emission lines (Sect. lITSl . All ions with uncertain column 
densities are indicated by * (see Sect. 13. 3> . 



ion 


log TVion 

o = 140 


log iVion 

a = 70 


V 


log? 


logf/ 


ion 


log iVion 

(7= 140 


log iVion 

a = 70 


V 


logC 


logU 


Civ* 


21.3 ± 0.1 


20.6 ± 0.4 


500 


-0.65 


-2.25 


ArXIII* 


19.6 ±0.6 


19.5 ± 0.5 


700 


1.5 


-0.10 


Cv 


20.6 ± 0.2 


21.2 ± 0.1 


530 


0.1 


-1.50 


ArXIV* 


< 20.0 


< 19.5 


700 


1.75 


0.15 


Cvi 


21.2 ± 0.1 


21.5 ± 0.1 


480 


1.0 


-0.60 


ArXV* 


< 20.0 


<20.1 


700 


1.95 


0.35 


Nv 


< 19.8 


< 19.8 


- 


-0.15 


-1.75 


ArXVI* 


< 21.7 


22.1 ±0.1 


800 


2.30 


0.70 


Nvi 


20.3 ± 0.2 


21.2 ±0.2 


550 


0.5 


-1.10 


ArXVII* 


<20.0 


<20.2 


800 


2.60 


1.00 


Nvil 


21.2 ± 0.2 


21.7 ± 0.1 


320 


1.35 


-0.25 


ArXVIII* 


< 21.2 


<20.7 


800 


3.05 


1.45 


Oiii* 


< 20.5 


20.2 ± 0.4 


400 


-1.7 


-3.30 


Caxi* 


< 18.8 


19.2 ± 0.4 


600 


0.70 


-0.90 


Oiv* 


20.3 ± 0.3 


20.2 ± 0.3 


400 


-0.85 


-2.45 


Caxil* 


< 19.4 
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600 


1.00 


-0.60 


Ov 


20.6 ± 0.2 
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-0.2 


-1.80 


CaXIII* 


< 19.6 
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700 


1.40 


-0.20 


Ovi 
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0.3 


-1.30 


CaXIV* 


19.7 ±0.3 
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1.60 


0.00 


Ovii 


21.83 ± 0.05 


21.94 ± 0.03 
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0.95 


-0.65 


CaXV* 


< 18.6 


< 19.2 
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2.00 


0.40 


O VIII 


22.31 ±0.04 


22.37 ± 0.03 
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1.65 


0.05 


CaXVI* 


<20.2 


< 19.8 


800 


2.25 


0.65 


Nevii* 


21.6 ± 0.1 


< 20.7 
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0.35 


-1.25 


CaXVII* 


20.0 ± 0.5 


< 19.8 
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2.40 


0.80 


Ne VIII* 


21.4 ± 0.1 


20.3 ± 0.8 


500 


0.9 


-0.70 


CaXVIII* 


20.2 ± 0.4 


< 20.2 


800 


2.60 


1.00 


Ne IX 


21.3 ± 0.1 


20.9 ± 0.2 


660 


1.5 


-0.10 


CaXIX* 


< 20.9 


< 20.5 


800 


2.80 


1.20 


Nex 


21.5 ± 0.2 


21.3 ±0.1 
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2.05 


0.45 


CaXX* 


< 20.9 


< 20.4 


800 


3.25 


1.65 
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21.4 ± 0.1 


< 20.3 
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0.25 


-1.35 
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< 23 
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-3.55 
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-3.20 
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-0.35 
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700 


1.7 


0.10 
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-1.00 


-2.60 


NaXI* 
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< 23 


700 


2.3 


0.70 


FeVII* 


19.9 ±0.6 


< 19.8 


500 


-0.50 


-2.10 


Mg VI* 


< 20.1 


< 19 


500 


-0.3 


-1.90 


Fe VIII* 


20.0 ± 0.3 


< 19.8 


500 


-0.10 


-1.70 


Mg VII* 


20.6 ± 0.9 


< 19.7 


500 


0.1 


-1.50 


FelX 


20.2 ± 0.4 


20.97 ± 0.08 


500 


0.20 


-1.40 


Mg VIII* 


<20.2 


< 19.8 


500 


0.5 


-1.10 


Fex 


20.4 ±0.1 


20.2 ± 0.1 


500 


0.40 


-1.20 


MglX 


20.3 ± 0.2 


20.3 ± 0.3 


560 


1.0 


-0.60 


FeXI 


20.1 ±0.2 


19.8 ± 0.2 


500 


0.60 


-1.00 


MgX* 


< 20.2 


< 19.4 


700 


1.4 


-0.20 


FeXII* 


19.7 ±0.5 


< 19.7 


600 


0.80 


-0.80 


MgXI 


20.9 ± 0.2 


20.4 ± 0.2 


500 


1.9 


0.30 


FeXIII 


19.5 ±0.6 


19.4 ± 0.4 


600 


0.95 


-0.65 


MgXII 


21.0 ±0.1 


20.8 ± 0.1 
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2.35 


0.75 


FeXIV* 


< 19.7 


19.3 ± 0.9 


700 


1.15 


-0.45 


Al VII* 


< 20.8 


< 19.9 


400 


0.05 


-1.55 


FeXV* 


< 19.6 


< 19.0 


700 


1.40 


-0.20 


Al VIII* 


21.1 ±0.2 


20.7 ± 0.4 


400 


0.45 


-1.15 


FeXVI* 


< 18.9 


< 18.6 


700 


1.40 


-0.20 


AllX* 


< 19.6 


< 19.3 


600 


0.9 


-0.70 


Fe XVII 


20.7 ± 0.2 


20.5 ± 0.2 


740 


1.80 


0.20 


AlXI* 


< 19.3 


< 19.5 


700 


1.7 


0.10 


Fe XVIII 


20.3 ± 0.2 


20.1 ± 0.2 


700 


2.05 


0.45 


AlXII* 


< 20.0 


< 19.9 


700 


2.05 


0.45 


FeXIX 


20.7 ±0.1 


20.5 ±0.1 


800 


2.25 


0.65 


AlxiII* 


< 20.3 


< 20.7 


800 


2.45 


0.85 


FeXX 


20.4 ± 0.2 


20.4 ± 0.1 


800 


2.40 


0.80 


Si VII* 


19.7 ± 1.0 


< 19.8 


400 


0.05 


-1.55 


FeXXI 


20.7 ± 0.2 


20.5 ± 0.2 


800 


2.60 


1.00 


Si VIII* 


19.7 ±0.3 


20.9 ± 0.1 


620 


0.4 


-1.20 


Fe XXII 


20.2 ± 0.2 


20.3 ± 0.3 


800 


2.75 


1.15 


Si IX* 


20.0 ± 0.4 


21.1 ± 0.1 


400 


0.75 


-0.85 


Fe XXIII 


20.5 ± 0.2 


20.4 ± 0.2 


800 


2.90 


1.30 


Six 


19.9 ±0.2 


20.0 ± 0.2 
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1.05 


-0.55 


Fe XXIV 


20.7 ± 0.2 


20.8 ± 0.2 


800 


3.10 


1.50 


Si XI 


19.7 ± 0.2 


19.8 ± 0.4 
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1.5 
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Nil* 


20.3 ± 0.2 
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< -4 
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Si XII 


20.6 ± 0.5 
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1.75 


0.15 
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0.60 
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0.6 


-1.00 
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20.3 ±0.1 
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400 


-0.95 


-2.55 


Sx* 


20.3 ± 0.3 
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600 


1.0 


-0.60 
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19.9 ±0.3 


20.1 ± 0.2 


500 


-0.50 


-2.10 


Sxi 


19.7 ±0.3 


19.7 ± 0.4 


560 


1.2 


-0.30 
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20.2 ±0.1 
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-0.05 


-1.65 
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620 


1.45 
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20.0 ± 0.3 


500 


0.40 


-1.20 
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20.1 ±0.4 
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600 
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-0.90 
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0.55 
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600 
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-0.70 
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2.4 
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NiXIII* 


20.0 ± 0.2 
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600 


1.10 


-0.50 


Sxvi* 


< 20.8 


21.0 ± 0.4 


800 


2.85 


1.25 


NiXIV* 


< 19.7 


19.5 ± 0.5 


700 


1.25 


-0.25 


ArIX* 


19.5 ±0.5 


19.6 ±0.3 


500 


0.4 


-1.20 


NiXV* 


< 18.9 


< 19.3 


700 


1.50 


-0.10 


Arx* 


19.9 ±0.3 


19.4 ± 0.3 


500 


0.8 


-0.75 


NiXVI* 


< 19.4 


19.7 ± 0.2 


700 


1.40 


-0.20 


ArXI* 


19.9 ±0.2 


19.7 ± 0.2 


600 


1.1 


-0.40 


Ni XVII* 


19.6 ± 0.4 


20.8 ± 0.3 


700 


1.60 


0.00 


ArXII* 


19.3 ±0.8 


19.2 ± 0.4 


600 


1.3 


-0.20 


Ni XVIII* 


19.7 ±0.3 


19.7 ± 0.3 


700 


1.70 


0.10 
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APPENDIX B 



In this Appendix we discuss those ions for which the column 
density may be affected by the effects described in Sect. 13.31 
Carbon, Oxygen and Neon: For C IV the strongest line at 41.42 
A coincides with the instrumental C I edge, which is deep 
and susceptible to calibration uncertainties. For O IV the wave- 
length of the dominant line in the X-ray range is not well known 
and the absorption line is weak. For Ne VIII the deepest line at 
88.10 A is at a rather noisy part of the spectrum, the other Une 
at 67.38 A is blended with an Al VIII Hne. 



Sodium and Aluminum: Both Na VII and Na VIII are only 
detected through lines above 60 A. The column density for 
Al VII is mainly determined from the 2pi/2 edge at 43.56 A and 
several lines above 60 A. A similar situation occurs for Al VIII, 
where the 2pi/2 and 2p3/2 edges at 51.36 A dominate in the 
column density determination. 



Silicon and Sulfur: The deepest lines for Si VIII and Si IX 
at 61.04 A and 55.30 A, respectively are saturated. For Si VIII 
this line is also blended. The saturation of these well detected 
lines explains the difference in column density obtained for the 
different velocity broadening assumed. For S X the strongest 
lines form a blend in the instrumental C I edge. 



Argon: The column density of Ar IX is determined from an 
absorption line which coincides with the instrumental C I edge 
and the C V forbidden line. Arx has only one detectable line 
at 37.45 A, which is at a rather noisy part of the spectrum. For 
Arxi the deepest line is blended with the —1040 km s^^ com- 
ponent of S X at 34.29 A. Ar Xll has only one weak line in the 
spectrum at 31.37 A, and for Arxill we detect only two weak 
Unes at 29.27 A and 27.44 A. The deepest ArXV absorption 
line is blended with the N Vll Lya line. The blending is further 
complicated by the possible substructure in the N VII Lya line 
due to the different outflow velocities. ArXVl is detected from 
only two edges in the spectrum at 15 and 25 A. 



Nickel: The detection of all nickel ions is uncertain. Nickel 
only produces weak absorption lines, most of them are blended 
with stronger iron absorption lines. Small errors in wavelength 
for these iron lines could be compensated in our model by ab- 
sorption from nickel. For Ni I the column density is determined 
from the observed short wavelength tail of the edge at 105.97 
A and a smafler edge at 14.42 A. For Ni II the 3s edge at 87.93 
A dominates the column density determination. The strongest 
Ni VI and Ni VII lines are blended by Fe XVI at 14.30 A and 
Fe XVIII at 14.26 A, respectively. The strongest Ni IX line is 
blended by an FeXX line at 13.73 A. The column density for 
Nixi was determined mainly from the 3p-edge at 38.62 A, 
close to the instrumental C I edge, and in a noisy part of the 
spectrum. The two strongest lines of Ni XII are blended with 
the Ne IX resonance line at 13.45 A. Ni XIII causes a slight de- 
pression in the spectrum due to many weak absorption lines. 
However, the detection is uncertain as the continuum spectrum 
of the MEG around the relevant wavelength of 13.27 A differs 
by about 5 % from the LETGS spectrum, a similar level as the 
expected depression. The deepest line for Ni XVII is blended 
with an FeXX at 12.90 A; Nixvill has several weaker lines 
that are blended with Fe XX at 12.60 A. 

Due to the uncertainties in the determined column densities 
for the above listed ions, we decided not to use them in the 
further analysis. 

APPENDIX C 

Here we present the LETGS, MEG and HEG spectra of 
NGC 5548. The dotted line indicate the zero flux level for the 
MEG and HEG spectra. The HEG spectrum was rebinned by a 
factor of 2 for plotting purposes, the MEG data was rebinned 
by a factor of 2 above 19 A. The Fe Ka line appears slightly 
shifted in wavelength between the HEG and the MEG spectra, 
due to calibration uncertainties or noise in the MEG spectrum. 
Between 4 A and 6.2 A there is a calibration mismatch between 
the HEG and the MEG spectra. The thick line is for the model 
with broadened emission lines, the thin line is the same model 
but without broadened emission lines. 



Calcium: Caxiv has only one unblended line at 21.13 A. 
However, the fit is driven by only one data point of the LETGS 
spectrum, making the column density uncertain. For all the 
other calcium ions the detection is very sensitive to the velocity 
broadening. If we leave a free, we mostly find upper limits to 
the velocity broadening of 30 km s^^, implying highly satu- 
rated lines and therefore unrealistically high column densities. 
The lower ionized states of calcium (Ca I to Ca Vlll) have edges 
in the spectra, which dominate the fit. 
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Fig. 12. Detail of the HEG, MEG and LETGS spectra. The dotted hnes indicate the zero flux level for the HEG and MEG spectra. 




Fig. 12. Continued. 



